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ABSTRACT 

We describe Sunyaev-Zel'dovich (SZ) effect measurements and analysis of the intracluster medium 
(ICM) pressure profiles of a set of 45 massive galaxy clusters imaged using Bolocam at the Caltcch 
Submillimeter Observatory. These clusters have a median mass of M500 = 9x 10 14 M and a median 
rcdshift of 0.42. We have used masses determined from Chandra X-ray observations to scale each 
cluster's profile by the overdensity radius -R500 and the mass-and-rcdshift-dcpcndcnt normalization 
factor -P5oo- We deproject the average pressure profile of our sample into 13 logarithmically spaced 
radial bins between 0.07i?5oo and 3.5i?5oo- We find that a generalized Navarro, Frenk, and White 
(gNFW) profile describes our data with sufficient goodness-of-fit and best-fit parameters (C500, a, j3, 
7, Pq = 1.18, 0.86, 3.67, 0.67, 4.29). We also use the X-ray data to define cool-core and disturbed 
subsamplcs of clusters, and we constrain the average pressure profiles of each of these subsamples. 
We find that given the precision of our data the average pressure profiles of disturbed and cool-core 
clusters are consistent with one another at R > 0.15i?5oo, with cool-core systems showing indications 
of higher pressure at R < 0.15i?5oo- In addition, for the first time, we place simultaneous constraints 
on the mass scaling of cluster pressure profiles, their ensemble mean profile, and their radius-dependent 
intrinsic scatter between O.I-R500 and 2.0i?5oo- The scatter among profiles is minimized at radii between 
~ 0.2i?5oo and ~ 0.5i?5oo, with a value of ~ 20%. The best-fit mass scaling has a power-law slope of 
0.49, which is shallower than the nominal prediction of 2/3 from self-similar hydrostatic equilibrium 
models. These results for the intrinsic scatter and mass scaling are largely consistent with previous 
analyses, most of which have relied heavily on X-ray derived pressures of clusters at significantly lower 
masses and redshifts compared to our sample. Therefore, our data provide further evidence that cluster 
pressure profiles are largely universal with scatter of ~ 20-40% about the universal profile over a wide 
range of masses and redshifts. Furthermore, and in agreement with recent results, our data show the 
consistency of X-ray and SZ measurements of the ICM given current observational precision. 
Subject headings: galaxies: clusters: general — galaxies: clusters: intracluster medium 



1. INTRODUCTION 

Massive galaxy clusters, the largest virialized sys- 
tems in the universe, appear to be exceptionally reg- 
ular objects. This is especially true at intermediate 
radii outside of the core (where complicated baryonic 
physics plays a large role), and inside the actively ac- 
creting outer regions (where non-equilibrium effects be- 
come significant). In these intermediate regions sim- 
ple, self-similar scalings based on hydrostatic equilib- 
rium and gravitational physics describe observations and 
simul ations quite well (jKaisedli98 6: Krav tsov fc Borganil 
[20121 1. Specifically, a range of observational results 
show approximately universal behavior among mass 
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and intracluster medium (ICM) profiles after scal- 
ing by characteristic overdensity radii such as i?5oc£3 



malizations (IVikh 
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For example, the cluster-to-cluster dispersion (intrin- 
sic scatter) observed in these scaled profiles out- 
side of the core regions and ins ide R^pg i s gen- 
erally ~ 10- 40% f or entropy (IPratt , et all [20061: 



Mm 
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perature (fvikhlinin et all I2006at IPratt et alt ' 120071: 
Leccardi fc Molendill2008i r and pressure fAlO. ISun et~aLl 
201 If ). In particular, both simulations and observa- 



tions indicate low cluster-to-cluste r dispersion in pres- 
sure profiles at in t ermediate radii ([Borgani et al.l 120041 : 
iNagai et all [20071 iPiffaretti fc Valdarninil 120081 A10, 

10 ^500 denotes the radius where the average enclosed mass den- 
sity is 500 times the critical density. Throughout this manuscript 
we define quantities at, or enclosed within, this characteristic ra- 
dius. 
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iPlagge et al1l2010t iBonamente et all 120121 P12). 

Historically, observational studies of the ICM pressure 
have relied almost exclusively on X-ray data. These 
data have provided precise constraints on the pres- 
sure profiles in the inner regions of clusters (R < 
i?50o), but the density-squared dependence of the X- 
ray surface brightness makes it difficult to study the 
ICM at large radii with X-rays. Although several X- 
ray results have extended beyond -Rroq in individual 
clusters or small sets of clusters (IGeorge et al.l [20091: 
Bautz et alj|2009t IReiprich et al.l 2009: Simio nescu et al.l 



20111: IWalker et alj|2012f >. with current X -ray instrumen- 
tation it is infeasible to extend such studies to large 
samples of clusters. Consequently, two groups have 
used a hybrid approach with X-ray data at small radii 
(R < -R500) and simulations at large radii (R > R500) 
in order to constrain t he average pressur e profile at all 
relevant radial scales ([Nagai et al.l 120071 A10). These 
X-ray and simulation-based results imply that cluster 
pressure profiles are approximately universal over a wide 
range of masses and radial scales, with low intrinsic scat- 
ter that is minim i zed near 5 R50 at < 20% (e.g., A10, 
INagai et al .1120071 1 Sun et al.ll201lh . 

The Sunyaev-Zerdovich (SZ) effect signal 
(|Sunvaev and Zel'dovichllT972D . which is proportional to 
the density of the ICM and therefore falls more slowly 
with radius compared to the X-ray brightness, can be 
exploited to study the ICM pressure at large radii. 
Although some initial studies using WMAP SZ data 
showed large inconsistencies with the establishe d X-ray 
results (|Lieu et al.ll2"00l iBielbv fc Shanks! [2001 . recent 
results have shown that SZ data from WMAP, Planck, 
and ground-based receivers provide a picture of the ICM 
that is consistent with X-ray measuremen t s at current 
observational precision (IPlagge et alj I2010t IMelin et all 
20111: \Planck early results Xll2011t IKomatsu et al.H2011t 
Bonamente et a.1.1 1201 2f ) . In particular, the South Pole 

to measure the average 
2i?5oo for a sample of 
finding results that 



Telescope (SPT) was able 
SZ pressure profile out to : 
15 clusters (|Plagge et all 12010(1 . 
were similar to previous X-ray /simulation results at 
small/large radii (e.g., the sample of 31 REXCESS 
clusters studied by A10). Recently, a combination of 
XMM-Newton X-ray data at small radii and Planck 
SZ effect data at large radii was used to constrain the 
average pressure profile out to 3i?5oo for a sample of 62 
Planck-selected clusters (P12). These results were again 
largely consistent with previous analyses, and the X-ray 
and SZ data agreed quite well in the overlapping region 
at intermediate radii. Altogether, X-ray and SZ data, 
along with simulations, are converging to a uniform 
picture of the average cluster pressure profile over a 
wide range of angular scales. 

This manuscript is arranged as follows. In Section[5]we 
describe our sample of 45 massive clusters, and in Sec- 
tion [3] we provide the details of our SZ and X-ray data 
reduction. We then present our method for deprojecting 
pressure profiles from our spatially-filtered SZ images in 
Section In Section [5J we describe parametric fits to 
these deprojected profiles, and compare the results of 
our fits to the results from a range of previous analy- 
ses. Then, in Section [51 we use a Gaussian process for- 
malism to simultaneously constrain the pressure-profile 
mass scaling, the ensemble mean profile, and the radius- 



dependent intrinsic scatter about this mean profile. Fi- 
nally, we provide a summary of our results in Section [7] 

2. CLUSTER SAMPLE 

Between 2006 November and 2012 March we used Bolo- 
cam to image the SZ signals from a sample of 45 clus- 
ters that have Chandra X-ray exposures (hereafter the 
Bolocam X-ray/SZ or BoXSZ sample). The Bolocam SZ 
observations and some general properties of the BoXSZ 
sample are summarized in Table [1] and the details of 
the Chandra X-ray data and analysis are given in Ta- 
bles [H and [3] The BoXSZ sample contains two previ- 
ously defined subsamplcs, the 25 object Cluster lensing 
and Supernova Surve y with Hubble (C1ASH) sample 
(iPostman et al.l l2012h and th e 12 cluster MACS high- 
z sample ( Ebeling et alj|2007[ l. The remaining clusters 
were selected in an ad hoc manner, with a general em- 
phasis on massive and/or high redshift systems. The 
clusters span the redshift range 0.15 < z < 0.89, with a 
median redshift of 0.42 and more than 60% of the sam- 
ple lying between 0.35 < z < 0.59. The clusters are 
among the most massive known, with a median X-ray 
derived mass of M500 = 9 x 10 14 M Q assuming a con- 
stant gas mass fraction (we computed these masses us- 
ing a reference h = 0.7, f2 m = 0.3 flat ACDM cosmology, 
and this same cosmology was used to calculate all other 
physical quantities presented in this manuscript). The 
clusters in the BoXSZ sample span a range of dynami- 
cal states, from rela xed systems with well defined cool- 
cores like Abell 1835 ([Peterson et al .11200 lUSchmidt et all 
I2001H to clusters undergoing major merger even t s like 
MACS J0717.5 (lEbeling et all I200H lEdge et al.l [200l 
iMroczkowski et al.ll2012j) . 

Based on previous results that show the projected X- 
ray lumi nosity ratio i s an accurate indicato r of cool-core 
clusters (|Mantzll2009llBohringer et al.ll2010f ). we define a 
cool-core subsample of the BoXSZ sample as those clus- 
ters with a projected X-ray luminosity ratio 



L(R < OMR. 



500 ) 



L(R < R, 



> 0.17. 



.5(1(1 ) 



(1) 



According to this definition, 17/45 of the BoXSZ clusters 
are cool-core systems (see Table [3]) . We note that these 
cool-core systems are in general at the low redshift end of 
our sample, with a median redshift of z = 0.36. We spec- 
ulate that this result is due at least in part to the fact that 
the cool-core fraction of clusters drops with increasin g 
redshift (e.g.. IVikhlinin et al|[2006al iSantos et al.ll2010jl . 
although selection effects may also play a role due to the 
ad hoc manner in which the BoXSZ sample was chosen. 

We define a disturbed subsample of the BoXSZ sample 
as those clusters with an X-ray ccntroid shift parame- 
ter of U7500 > 0.01 (see Section [3721 for a full description 
of how we compute u>5oo)- The centroid shift param- 
eter is widely used to classify disturbed systems (e.g. 



Maughan etail l2008t IPratt et all 120091: IMaughan et all 



tug. 

2012), and we ado pt the same threshold (u>5oo > 0.01) as 



Pratt etTaTI (|2009( 1. Based on this criteria, 16/45 clusters 
in the BoXSZ sample are disturbed, and the disturbed 
systems are generally at the high redshift end of the full 
sample with a median redshift of z = 0.52. This result is 
not surprising given that the cool-core systems are gen- 
erally at low redshift and cool cores are a good indicator 
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TABLE 1 

Bolocam X-ray/SZ (BoXSZ) Cluster Sample 



cluster redshift RA dec obs. time noise peak S/N 
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Note. — The Bolocam X-ray/SZ (BoXSZ) cluster sample of 45 objects. The columns give 
the name, redshift, X-ray centroid coordinates (J2000), total Bolocam integratio n time, RMS 
noise l evel of the SZ images, and peak S/N in the optimally filtered images (sec Savers et al. 
2012b for details of how the peak S/N and op timal filter are dete rmi ned). The superscr ipts 
denote the reference fo r th e redshifts, with a) | Mantz et al.l J2010I), b ) lAllen et all 1120081), c) 
Maughan eFaLI (120121), d) IRichard et al.l 1 12010) , e) Uee et all H2007n , fT TChristensen et alj 
H2012I ). and o) ITran etUl l|2007l ) . 



that the cluster is not disturbed (only 2 /45 clusters in the 
BoXSZ sample have both a cool core and are disturbed). 

This redshift asymmetry between the cool-core and dis- 
turbed subsamples within the BoXSZ sample has an im- 
pact on our analysis of the average pressure profile of the 
sample. As we describe in Section 0] we scale the radial 
coordinate of each cluster by -R500, since the shape of 
the pressure profiles is expected to be self-similar after 
scaling by this radius. However, the angular dynamic 
range of our data is limited by Bolocam's 58 arcsec full- 
width at half-maximum (FWHM) point-spread function 
(PSF) and by the 14 arcmin size of the Bolocam images 
(i.e., our SZ data are sensitive to a minimum radius of 
~ 29 arcsec and a maximum radius of ~ 10 arcmin). 



Since the value of -R500 in physical units is fairly con- 
stant over our sample (~ 1.5 Mpc), this means that the 
value of i?5oo varies significantly in angular size over the 
redshift range of our sample, from 2 arcmin to 9 arcmin. 
This can be clearly seen in Figure [U where we show the 
Bolocam integration time as a function of scaled radius 
for the BoXSZ sample. On average, we obtained longer 
integrations for the preferentially high redshift disturbed 
clusters, and this compensates for their smaller angu- 
lar size compared to the preferentially low redshift cool- 
core clusters to produce nearly equal integration times 
at R < i?5oo- However, the increased integration time 
and smaller angular sizes of the disturbed clusters re- 
sults in significantly more integration time outside -R500 
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TABLE 2 

Details of the Chandra ACIS observations used in this 
work. 



cluster date ObsID mode exp. (ks) 
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compared to the cool-core clusters. Consequently, the 
average SZ pressure profile of the BoXSZ sample is con- 
strained by an above average amount of data from dis- 
turbed systems at larger radii. At these large radii, data 
and simulations indicate that there is little or no dif- 
ference in the pressure profiles based on morphological 
classification, so we ex pect that any resulting bias in our 
resul t s will be minor dBorgani et al.1 120041 : iNagai et al.1 
[20071 IPiffaretti fc Valdarninil I2008T A10'). We address 
this issue in more detail in Section [5j where we exam- 
ine the average pressure profiles of our cool-core and 
disturbed subsamples separately and find that they are 
consistent outside of ~ 0.15i?5oo given our measurement 
uncertainties. 

3. DATA REDUCTION 
3.1. Bolocam 

Bolocam is a 144-clement bolometric imaging pho- 
tometer, and from 2003-2012 it served as the long- 
wavelength facility camera for the Caltech Submillimctcr 
Observatory (CSO). Bolocam covers an eight-arcminutc- 
diameter circular fi eld of view (FOV) and has PSFs wit h 
58 arcsec FWHMs ([Glenn et al.lll998t IHaig etllll200l . 
The SZ-emission- weighted band center of our data is 
140 GHz. All of the cluster images were obtained by 
scann ing the CSO in a Lissajous pattern (jKovacs et al.l 
120061 ) , with an amplitude of 4 arcmin and an average scan 
speed of ~ 4 arcmin/sec. These scans result in images 
with tapered coverage extending to a radius of ~ 12 ar- 
cmin, with the coverage dropping to half its peak value 
at a radius of ~ 5 arcmin. For ease of analysis, we have 
made 14 x 14 arcmin square maps for each cluster. 




0.1 1.0 

radius (R 500 ) 

Fig. 1. — Total Bolocam integration time for BoXSZ clusters as 
a function of scaled radius. The black line denotes the full BoXSZ 
sample, the red line denotes the disturbed subsample, and the blue 
line denotes the cool-core subsample. 

Our Bolocam data reduction largely followed th e pro- 
cedure described in detail in I Savers et al.l (| 2 1 If ) . and 
we therefore briefly summarize that procedure below. 
First, we use frequent observations of bright compact ob- 
jects to obtain pointing corrections accurate to 5 arcsec. 
Additionally, we made nightly observations of Uranus, 
Neptune, and/or other seconda ry calibrators to obtain 
flux calibration accurate to 5% ( Griffin and Ortorj||1993t 
ISande ll 1994; Save rs et al.H2012af ). To remove noise from 
atmospheric fluctuations, we subtract the FOV-average 
signal at each time sample in the time-ordered data 
(TOD) and also high-pass filter the TOD at a character- 
istic frequency of 250 mHz. This process also attenuates 
the astronomical signals in our data, and we characterize 
this filtering as follows. First, for each observation, we 
insert a model cluster profile in our TOD, process these 
modcl-plus-data TOD through our reduction pipeline, 
and create a map. We then subtract the data-only map 
to produce a noiseless image of the model after going 
through our data processing pipeline. The result is com- 
pared to the original input model to obtain a complex- 
valued two-dimensional map-space transfer function. 

To characterize the non-astronomical noise in our im- 
ages we form jackknife realizations of the data by multi- 
plying a randomly selected subset of half of the data by 
— 1 prior to binning the data into a map. For each cluster 
we formed 1000 such jackknife maps. To each of these 
maps we then added a Gaussian random realization of 
the 140 GHz astronomical sky based on the power spec- 
trum measurements made by the SPT, w hich cover all 
of the angular scales probe d by our data (jKeisler et al.l 
I20T1 IReichardt et al.l 120121) . Each of these 1000 astro- 
nomical signal realizations was processed through our 
data reduction pipeline so that it was filtered identi- 
cally to our real data. We have verified that the above 
noise model is statistically equiva lent to measureme nts 
of blank sky made with Bolocam (ISavers et al.ll20fTI ). 

In addition, for 11 of our clusters we have subtracted 
individual bright point sources selected from the NVSS 
1.4 GHz catalog dCondon et all 119981). We refer the 
reader to ISavers et al.l (|2012cf) for a full description of 
these sources, most of which are near the cluster cen- 
ters. Since all of these sources are below our 140 GHz 
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± 
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1 
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CL J1226.9 
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.00 
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.8 


± 


1. 
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1 


.3 


11 


.84 


± 


3 


.21 





.10 


± 





.02 





.95 


± 





.31 







Note. — Relevant X-r ay derived propertie s of the BoXSZ cluster sample. The first four columns provide i?500i ^500 1 -^500 1 an d fcT, 
computed as described in Mantz et al. ( 2010T). The fifth column gives P50O1 computed from M500 and 2: according to Equation[2] given 
in Section [4] The sixth column gives the ratio of the projected X-ray luminosity within 0.05i?5oo to the X-ray luminosity within -R500 
(Mantz 2009). The seventh column gives the centroid shift parameter within K50O1 computed according to the procedure described in 
Maughan ct al. (2008, 2012). Wc denote cool -core clusters as those having an L r at > 0.17, and we denote disturbed clusters as those 
having W500 > 0.01. 



detection limit, we have extrapolated spectral fits to 1.4 
and 30 GHz data. In almost all cases, the uncertainty 
on these extrapolated flux densities is ~ 30%, limited 
by the intrinsic scatter in the extrapolation. We sub- 
tracted all of the sources with extrapolated flux densities 
> 0.5 mJy. This source brightness threshold was chosen 
to ensure that contamination of the cluster signal from 
point sources is < 1%. To account for our uncertainty 
in the flux density of these subtracted point sources, we 
add a model of each point source, multiplied by a random 
value drawn from a Gaussian distribution described by 
our uncertainty on the extrapolated source flux density, 
to each of the 1000 noise realizations for a given cluster 
observation. 

Furthermore, we detect a total of 6 bright point sources 
in our 140 GHz data, all of which were subtracted for 



this analysis (|Savers et al.H2012cT) . For these sources, we 
refit the point source model to each of our 1000 jackknifc 
realizations and added a point source template to each 
of our 1000 noise realizations based on the dispersion of 
these fits. 

3.2. X-ray data 

X-ray luminosities, temper atures and masses for the 
BoXSZ clusters appearing in iMantz et all ()2010h (here- 
after M10) arc taken from that work. For the other clus- 
ters, these quantities were derived from archival Chandra 
data following the same procedure as was used in M10, 
and we refer the reader there for full details. Briefly, 
the archival data were reprocessed using ciac0 (version 

11 http://cxc.harvard.edu/ciao/ 
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4.1.1; CALDB 4.1.2), including removal of bad pixels, 
corrections for cosmic ray afterglows and charge transfer 
inefficiency, and application of standard grade and status 
niters, using appropriate time-dependent gain and cali- 
bration products. Soft-band surface brightness profiles 
were extracted and were scaled by a global factor to agree 
with the final ROSAT flux calibration. These profiles 
were then used to derive cluster luminosity, projected lu- 
minosity ratio (Equation [J), and gas mass. Values of 
total mass and -R500 were derived for a reference h = 0.7, 
TO = 0.3, flat ACDM cosmology using the derived gas 
mass pr ofiles and the univ ersal gas mass fraction mea- 
sured bv lAUen et all ((2008V Spectra were extracted from 
an annulus covering radii between 0.15i?5oo and -R500 and 
fit in xspecQ to provide a global temperature measure- 
ment. For the lowest redshift cluster (Abell 2204), the 
luminosity and gas mass analyses used ROSAT Position 
Sensitive Proportional Counter data, since the Chandra 
field of view does not comfortably encompass -R500, and 
the average temper ature measured from ASCA data by 
iHorner et all (fl999t) was adopted. 

In addition, we used the Chandra data to constrain 
the centr oid shift parameter w^n p. This analysis was de- 
scribed in lMaughan et all (|2012f ). and we refer the reader 
to that manuscript for additional details. Centroid shifts 
measure the standard deviation of the projected sepa- 
ration between the X-ray peak and the centroid as a 
function of proje cted radius R n < -Rbdo- Following the 
same approach as lPoole et al.1 ()2006l ). w^oo was computed 
from a series of circular apertures with initial and final 
radii R p = R500 an d R P = 0.05i?,5oo respectively, de- 
creasing in size by O.O5A500 in each iteration. To this 
end, we used background-subtracted images, appropri- 
ately divided by the exposure map to eliminate instru- 
mental artifacts such as chip gaps and vignetting. All 
point sources were excluded from the analysis; however 
extended sources were left untouched as these may be 
associated with some of the cluster substructure. 

4. PRESSURE DEPROJECTIONS 

In order to determine physical pressures from our SZ 
data (which are measured in units of a CMB fluctu- 
ation temperature) we use th e equations described by 
iSunvaev and Zel'dovichl (11972T) wit h the relativistic cor- 
rections given bv lltoh et al For the relativistic 
corrections, we have assumed that the ICM of each clus- 
ter is isothermal, with a temperature equal to the spec- 
troscopic X-ray temperature given in Table [3l Although 
the true temperature profiles are expected to vary by 
factors of ~ 2 over the radii probed by our data, we 
note that the relativistic corrections are < 10% for these 
clusters in our observing band. Therefore, even if the 
ICM temperature varies by a factor of two over the radial 
range probed by our data, then our isothermal approxi- 
mation will produce less than a 5% bias on the pressures 
that we derive. 

Due to the high-pass filtering applied to the Bolocam 
data to remove atmospheric noise, it is not possible to di- 
rectly deproject pressure profiles from the standard Bolo- 
cam images. We have demonstrated the ability to de- 
convo lve the effects of this high pass filter dSavers et al.l 
120111 ). and we are able to produce unbiased images that 

12 http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/ 



could in principal be deprojected. However, to prevent 
unphysical numerical artifacts from appearing in the de- 
convolved images, we would need to reduce the size of 
our images from 14 to 10 arcmin. Since we are inter- 
ested in information about the pressure profiles at large 
radii, we do not want to accept this loss of information. 

Consequently, we compute our deprojected pressure 
profiles as follows. First, we select a set of discrete radial 
points (referred to as radial bins) at which we would like 
to compute the spherically symmetric deprojected pres- 
sure. Next, we create a smooth and continuous pressure 
profile by connecting these radial bins with a power-law 
interpolation. This smooth pressure profile is then pro- 
jected into a two dimensional SZ model. For the projec- 
tion, we assume that the profile is equal to zero outside 
of 5i?5oo, and we use a power- law extrapolation to esti- 
mate the pressure profile beyond our largest radial bin 
and inside 5i?soo- We then filter this SZ model with both 
the signal transfer function of our data and the Bolocam 
PSF, and then compare the filtered SZ model to our data 
map. 

We determine the best fit pressure deprojections us- 
ing the MPFIT FUN generalized l east squares (GLS) fit- 
ting algorithm (jMarkwardtl 120091 ) under the simplifying 
assumption that the noise covari ance matrix of our SZ 
maps is diagonal. As described in lSavers et al.l (|2011l) . a 
diagonal noise covariance matrix is a very good, but not 
perfect, description of our SZ map data. Consequently, 
we compute all of our uncertainties using the 1000 dif- 
ferent noise realizations described in Section 13.11 First, 
we add a cluster model, equal to the best fit deprojec- 
tion of our data, separately to each of these 1000 noise 
realizations. We then fit a pressure deprojection to each 
of these 1000 model-plus-noise realizations, and quantify 
all of our uncertainties based on the spread of these 1000 
fits. We have verified that our GLS algorithm adequately 
samples the parameter space by showing that it is robust 
to our choice of initial conditions, and that we are able 
to recover an input cluster candidate with minimal bias 
(see Figure [2]) . 

Throughout this work we compare and average the 
pressure profiles from multiple clusters. Therefore, prior 
to any deprojection, we have scaled the radial coordinate 
of each cluster by -R500 and the pressure amplitude by 

P„„= (3.68*10-^) (J^)" E(s r,( 2 ) 

where E 2 (z) = Um(^+z) 3 +^Ia and ap = 2/3 is the nom- 
inal scaling p redicted by self-sim ilar hyd rostatic equilib- 
rium models (|Nagai et al.ll2007l A1O)0 

We have determined the average pressure profile of the 
full BoXSZ sample (along with the disturbed and cool- 
core subsamplcs) via a simultaneous joint deprojection of 
multiple clusters. In this approach, a single deprojected 
profile in units of P500 and -R500 is constrained by the 
data from an arbitrary number of clusters using the GLS 
algorithm described above. Pressure deprojections for 

13 Unlike e.g. A10 or P12, our joint deprojections do not include 
any corrections to the value of ap , mainly because such corrections 
have negligible effects on our results. For example, correcting the 
value of ap from 2/3 to our best-fit ap = 0.49 found in Section [6] 
causes the pressure values in the individual joint deprojection bins 
for the full sample to change by < 2%. 
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Fig. 2. — The GLS fitting bias divided by the uncertainty for 
each radial bin in the joint pressure deprojection of the full BoXSZ 
sample. The black diamonds show the difference between the input 
pressure profile and the average recovered pressure profile using 
our GLS algorithm. The red triangles and green squares show the 
additional bias associated with choosing different starting values 
for the pressure bins relative to our default starting values equal to 
the best-fit model from A10 (For the red triangles all of the starting 
values were set equal to P500 , and for the green squares all of the 
starting values were set equal to a randomly drawn value between 
and lOPsoo)- In all cases, the bias is quasi- negligible when added 
in quadrature with the uncertainty (note that the typical S/N per 
bin is ~ 5, so the absolute bias is ~ 2%). This result shows that 
our fitting method is approximately unbiased and that it is robust 
to our choice of initial conditions (i.e., it adequately explores the 
parameter space). 

the full BoXSZ sample, the cool-core subsample, and the 
disturbed subsample are shown in Figure[3]and numerical 
values are given in Table |H For the deprojection of the 
full BoXSZ sample, wc have used 13 approximately log- 
arithmically spaced bins between 0.07i?5oo and 3.5i?5nn- 
Due to the smaller number of clusters in the two sub- 
samples, we have used 7 bins spanning the same radial 
range in each case. We note that, for many of the clus- 
ters in our sample, particularly the preferentially higher 
redshift disturbed clusters, the innermost radial bin(s) 
is(are) inside of Bolocam's PSF half-maximum radius of 
29 arcsec (in the most extreme case 0.07i?soo is — 8 arc- 
sec). However, we have shown that even for the disturbed 
subsample we are able to obtain unbiased fits to the in- 
nermost bins for deprojections of model input clusters, 
indicating that our data arc sensitive to the shape of the 
pressure profile at these sub-PSF-sized scales. 

We determine the covariance matrix of our dcprojcctcd 
pressure profiles directly from the set of 1000 deprojec- 
tions of modcl-plus-noisc data. The off-diagonal elements 
of the covariance matrix are in general non-zero and pos- 
itive at large radii, indicating significant positive corre- 
lations between most of those deprojection bins (see Fig- 
ure S]) . This can be understood as a consequence of the 
high-pass filtering applied to our data, which results in 
SZ images that are sensitive to the large-scale shape of 
the pressure profile (i.e., dP/dR) but not its absolute nor- 
malization. In contrast, the high-pass filtering does not 
have a significant effect on the deprojection bins at small 
radii, and adjacent bins at small radii tend to have strong 
negative correlations. These large correlations, which ex- 
ist for both the jointly and independently dcprojcctcd 
profiles, must be accounted for in any interpretation of 
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Fig. 3. — Pressure deprojection for the full BoXSZ sample (thick 
black line), the disturbed subsample (thick red line), and the cool- 
core subsample (thick blue line). Plotted as thin grey lines are the 
pressure deprojections for each of the 45 clusters in the BoXSZ 
sample. Each individual cluster was deprojected into 5 radial bins 
located at 0.5, 1.0, 2.0, 4.0, and 8.0 arcmin in order to fully sam- 
ple the spatial dynamic range of our images. One of the clusters is 
shown as a darker grey line with error bars at each deprojection bin 
to indicate the typical uncertainty on each of the individual cluster 
deprojections. A significant amount of the variation in the indi- 
vidual cluster profiles is due to measurement uncertainty, but we 
do find an additional clustcr-to-clustcr dispersion that is described 
in Section [6] Note that, in all cases, the correlations between bins 
are large and varied and must be accounted for in any fit or inter- 
pretation of the data. 
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Fig. 4. — Correlation matrix for the pressure deprojection of the 
full data set into 13 radial bins from 0.07-Rsoo to 3.5-Rsoo (see Ta- 
ble [4] for the exact radius of each bin). At large radii the adjacent 
bins have significant positive correlations due to the high pass filter- 
ing applied to the data, while at small radii the high-pass filtering 
has little effect on the data, and consequently there are significant 
anti-correlations between adjacent bins at those radii. 

our results. 

In addition to the joint deprojections described above, 
we also dcprojcctcd the pressure profiles of each cluster 
individually using an identical technique. However, due 
to the significantly varied range of scaled radii probed 
for each cluster, deprojecting all of the clusters at the 
same set of scaled radii is not ideal. Consequently, we 
have dcprojcctcd the individual clusters at fixed angu- 
lar radii to ensure that they adequately constrain the 
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TABLE 4 

Pressure Deprojections 



radius 


all clusters 


disturbed 


cool-core 


P500 


Psoo 




Psoo 



0.07 


(1.79 ± 0.30) 


X 


10+ 1 


(0.87 ±0.16) 


X 


10+ 1 


(2.07 ±0.14) 


x 10+ 1 


0.14 


(6.32 ± 1.48) 


X 


10° 












0.20 


(5.31 ± 1.11) 


X 


10° 


(4.81 ± 0.43) 


X 


10° 


(4.41 ± 0.29) 


x 10° 


0.28 


(2.68 ±0.38) 


X 


10° 












0.40 


(1.81 ±0.23) 


X 


10° 


(1.86 ±0.15) 


X 


10° 


(1.74 ±0.15) 


x 10° 


0.50 


(1.11 ±0.13) 


X 


10° 












0.63 


(6.41 ± 1.01) 


X 


10- 1 


(6.69 ±0.94) 


X 


10- 1 


(6.64 ± 0.93) 


x lO" 1 


0.80 


(4.82 ± 0.69) 


X 


10- 1 












1.00 


(2.44 ±0.61) 


X 


10- 1 


(3.22 ±0.62) 


X 


10- 1 


(2.63 ±0.56) 


x lO" 1 


1.26 


(1.87 ±0.45) 


X 


10- 1 












1.60 


(0.81 ±0.39) 


X 


10- 1 


(1.11 ±0.40) 


X 


10- 1 


(< 0.62) x 


10- 1 


2.00 


(4.99 ±3.72) 


X 


Kr 2 












3.50 


(3.60 ±2.63) 


X 


10~ 2 


(2.64 ±2.12) 


X 


10~ 2 


(< 8.13) x 


10~ 2 



Note. — Deprojected pressure profiles for the full BoXSZ sample, the dis- 
turbed subsamplc, and the cool-core subsample. The error bars give the square 
root of the diagonal elements of the covariance matrix. The best-fit pressure 
in the two outermost bins of the cool-corc dcprojection is consistent with zero, 
and so we instead list 68% confidence level upper limits on those values. 



profile over the scales probed by our data (0.5, 1.0, 2.0, 
4.0, and 8.0 arcmin). These individual deprojections are 
shown in Figure G2 and we describe our analysis of them 
in Section [51 

5. PARAMETERIZED PRESSURE PROFILES 

One method for characterizing the gross behavior of 
pressure profiles in our sample (or a subsample) is to fit 
a parameterized function to the data. We fit models to 
our deprojected data, rather than directly to the map 
data, because the former takes several minutes per fit 
compared to several days for the latter. As we justify 
below, the two fits are equivalent, and therefore there 
is no bias or loss of information involved in fitting the 
deprojected data rather than the map data. 

To perform the parametric fits to our deprojected data, 
we make the simplifying assumption that the noise is 
Gaussian and therefore fully described by its covariance 
matrix. We applied Mardia's test of multi variate norma l- 
ity to confirm this assumption is correct (Mardia 1970). 
Mardia's skewness statistic has a limiting distribution 
equal to the x 2 distribution, and we find a x 2 P er degree 
of freedom (DOF) equal to x? d = 465/455; Mardia's 
kurtosis statistic has a limiting distribution equal to the 
unit normal distribution, and we find a value of —0.96. 

As a final check, we directly fit a single generalized 
Navarro, Frenk, and White (NFW, iNavarro et"al~lll996| ) 
model to our map data, allowing all five parameters to 
vary. This model is described by 



P(X) 



Pn 



(C500WI + (C 50 o*) Q p-^/ Q ; 



where P{X) is the scaled pressure profile (in units of 
-P500), X = i?/i?5oo, Pq is the normalization, C500 sets 
the radial scale, and a, (3, and 7 describe t he power law 
slope at moderate, large, and small radii (jNagai et al.l 
|2007| ). Similarly to our deprojection procedure, we then 
determined our parameter uncertainties using fits to 1000 
separate noise realizations. For comparison, we per- 
formed a Markov chain montc carlo (MCMC) fit of the 
same model to our deprojected profile assuming a Gaus- 



sian covariance matrix, and we find that the best-fit pa- 
rameters are consistent within < 5% relative to their un- 
certainticsQ Therefore, we conclude that fitting para- 
metric profiles to our deprojected data is equivalent to 
fitting parametric profiles directly to our maps. Con- 
sequently, the six separate gNFW fit results in Table [5j 
along with the /3-model fit described below, were all com- 
puted using the deprojected data rather than the full 
map data, the latter of which would have taken several 
weeks of computing time compared to roughly an hour 
of computing time for the former. 

We find that our data do not constrain all five of the 
gNFW fit parameters within the physically allowed re- 
gion. Specifically, our best-fit outer slope of j3 = 2.03 im- 
plies an infinite total pressure^ As a result, we refit the 
data fixing the concentration parameter to C500 = 1-18 
(the best fit value found by A10), obtaining a physically 
allowed outer slope of j3 = 3.67. The quality of this 
four-parameter fit is good (x? e d = 1-0 for 9 DOF), indi- 
cating that this four-parameter gNFW fit describes our 
data with a sufficient goodness of fit. We note that the 
actual pressure profile over the radial range constrained 
by our data (0.07i?5oo < R < 3.5i?5on) for this four- 
parameter gNFW fit is only slightly different from the 
five-parameter fit, indicating that the unphysical outer 
slope in the five-parameter gNFW fit is likely due to the 
finite radial extent of our data. Therefore, the best-fit 
gNFW model of the BoXSZ sample is given by 

[C 500 ,a,P, 7 ,P ] = [1.18,0.86,3.67,0.67,4.29]. 

Our choice to fix the value of C500, rather than one of 
the other parameters, was motivated by the fact that 

14 We note that a Kolmogorov-Smirnov (KS) test shows that the 
exact noise distributions are inconsistent for 3/5 parameters, due 
to some extreme outlier values in the direct map fits that are not 
present in the MCMC fits to the deprojected data. However, due 
to the large degeneracies in the gNFW model, even these extreme 
parameter outliers in the fits to our map data produce pressure 
profiles consistent with the average profile over the range of radial 
scales probed by our data given our measurement uncertainties. 

15 Here, and throughout this manuscript, we present best-fit 
parameter values without error estimates due to the large degen- 
eracies between parameters, as illustrated in Figure [5] 
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TABLE 5 

gNFW Fits to the Bolocam Cluster Sample 



.Mm 



DOF PTE 



all clusters 

1.81 1.33 

1.18 1.05 

1.18 0.86* 

3.19* 3.28* 



4.13 
5.49 
3.67* 

2.03* 



disturbed clusters 
1.18 0.90* 5.22* 

cool-core clusters 
1.18 2.79* 3.51* 



0.31 
0.31 
0.67* 
1.10* 



0.02* 



1.37* 



6.54* 
7.82* 
4.29* 

3.07* 



17.28* 



0.65* 



1.7 
2.7 
1.0 
0.9 



3.1 



1.7 



12 
12 
9 

8 



0.00 
0.00 
0.44 
0.52 



0.02 



0.15 



Note. — gNFW fits to the deprojected profiles computed for 
the BoXSZ sample. The upper rows show fits to the deprojection 
of the full sample, and show parameter values and fit qualities 
when different numbers of parameters are varied. The lower two 
rows show the best-fit four-parameter models for the disturbed 
and cool-core subsamples. In all cases, stars denotes parameters 
that were varied in the fit. In the first two rows, we fit the models 
of P12 and A10, varying only the normalization. In all subsequent 
rows, the fixed parameters were set to the values found in A10. We 
find that varying four parameters provides a sufficient goodness 
of fit to the data, which is only marginally improved when all 
five parameters are varied. However, varying all five parameters 
results in a profile with an unphysically small outer slope (/3 = 
2.03), and we therefore take the four-parameter fit as the best 
description of our data. 
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Fig. 5. — Two-parameter confidence regions (68%) and one- 
parameter likelihoods for a gNFW fit to our joint deprojected pro- 
file of the full dataset. From left to right and top to bottom the 
plots show a, /3, 7, and Po with fixed C500 = 1-18. The large de- 
generacies between the fit parameters are clearly seen, along with 
the corresponding lack of constraining power on any individual pa- 
rameter. However, as shown in Figure|5] the overall pressure profile 
is tightly constrained. 

all of the possible four-parameter gNFW fits have sim- 
ilar fit qualities, but any combination that varies both 
C500 and j3 results in an unphysical outer slope. The 
set of two-dimensional confidence regions for the four- 
parameter gNFW fit is illustrated in Figure [5J and the 
strong parameter degeneracies mentioned above are evi- 
dent in these plots. 
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Fig. 6. — gNFW parameterized fits to the Bolocam data, varying 
four parameters of the gNFW model (C500 was fixed to 1.18). The 
bands indicate the maximum and minimum pressure values as a 
function of radius, bounding 68% of the MCMC fits. The plot 
shows fits to the full BoXSZ, the disturbed subsample, and the cool- 
core subsample. The pressure profile appears to be independent of 
cluster morphology at R > 0.15i?500i but the cool-core clusters 
have higher pressures than the disturbed clusters at smaller radii. 



We then fit the same four-parameter gNFW model 
to both the disturbed and cool-core subsamples of the 
BoXSZ cluster sample (see Table [5]). In good agreement 
with previous results from A10 and P12, we find consis- 
tent average pressure profiles, given our measurement un- 
certainties, between these subsamples (and the full sam- 
ple) at all intermediate and large radii (R > 0.15i?5oo)- 
At smaller radii, the profiles clearly diverge, and at 
those radii the cool-core clusters have a hi ghe r pressure 
than the disturbed clusters (see Figure [S]) E3 Although 
the overall pressure profiles for the full BoXSZ sample, 
the cool-core subsample, and the disturbed subsample 
are quite similar at most radii, we note that the best- 
fit gNFW parameters are quite different, further em- 
phasizing the large degeneracies between these param- 
eters. Furthermore, the fit quality of the gNFW model 
is slightly worse for both of the two subsamples compared 
to the full sample, with probabilities to exceed (PTEs) 
of 0.02 and 0.15 for the disturbed and cool-core fits com- 
pared to a PTE of 0.44 for the full-sample fit. 

Due to the large degeneracies between the parameters 
in the gNFW model, it is difficult to quantify the differ- 
ences between our best-fit gNFW model and those found 
in previous analyses via a dir ect co mparison of the fit 
parameters (jNagai et al.ll2007l A10. iPlagge et al.ll2010l 
P12). Consequently, we have compared the profiles re- 
sulting from our gNFW fits to the gNFW profiles found 
by A10 and P12 over the approximate radial range con- 
strained by all three datasets (0.05i?5oo < R < 4i? 50 o, see 
Figures [7] and [8]) . We in general find excellent agreement 
between our pressure profiles and those found in these 
previous analyses, regardless of morphological classifica- 

16 As we described in detail in Section [2] there is an asym- 
metry in the redshift distributions of the cool-core and disturbed 
subsamples. Consequently, we cannot definitively rule out redshift 
evolution as the cause of the discrepancy in the pressure profiles at 
small radii. However, given that A10 observed a similar difference 
in the inner pressure profiles of disturbed and cool-core systems 
for a sample of low-redshift clusters, the differences we observe are 
likely due to morphology rather than redshift evolution. 
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Fig. 7. — gNFW parameterized fits to the BoXSZ sample, varying four parameters of the gNFW model. From left to right the three plots 
show the Bolocam fit to the full sample, disturbed subsample, and cool-core subsample as points with error bars and as shaded regions, 
with the best-fit parameterizations given in A10 and P12 overlaid as thin and dashed lines (P12 did not fit a disturbed subsample, so we 
overlay their non-cool-core fit in the center plot). The A10 fits relied on the REXCESS sample of 33 low- z clusters (z < 0.2) observed with 
XMM-Newton within -R500 and results from simulations outside -Rsoo- The P12 fits relied on a sample of 62 Planck selected clusters at 
(z) ~ 0.15, and used XMM-Newton data to constrain the inner portion of the profile and Planck data to constrain the outer portion of the 
profile. Our fits use Bolocam SZ data for a sample of 45 higher redshift clusters (0.15 < z < 0.89). 
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Fig. 8. — The ratio of our best-fit four-parameter gNFW fits to the best-fit four-parameter gNFW fits from A10 (left) and P12 (right). 
In both cases the agreement is generally good in the regions that are well constrained by all three datasets (O.I-R500 5; R $ l-O-Rsoo). 
However, the fit to the full BoXSZ sample shows hints of higher pressure than the A10 fit at both large and small radii, and hints of higher 
pressure than the P12 fit at small radius. 

tiorTH (e.g., the cool-core profile from our analysis is in 
good agreement with the cool-core profile of A10). We do 
note that our disturbed and cool-core systems indicate 
slight differences within 0.3i?soo compared to the corre- 
sponding results of A10, and our overall average profile 
indicates slightly higher pressures at R < O.I-R500 and at 
R > I.O-R500 compared to the results of A10. Our over- 
all average profile also shows higher pressure at small 
radii compared to the results of P12. Our results are 
therefore more similar to simulation-derived results at 
small radii, which also s how higher pressures than found 
by A10 and P12 (A10, iBorgani et al] [200l iNagai et al.l 
I2007I : IPiffaretti fc Valdarninill2008l P12). 

The overall good agreement between our best-fit 
gNFW profiles and the results derived in previous anal- 
yses provides further evidence that the average clus- 
ter pressure profile is approximately universal (at least 
within our measurement uncertainties on the average 
profile, which are ~ 10-20% inside i?5oo)- This is es- 
pecially true given the large differences in the median 

17 P12 only present results for cool-core and non-cool-core sub- 
samples, and we therefore take their non-cool-core results to be 
representative of disturbed systems. 



rcdshifts ((z) = 0.12, 0.15, and 0.42), median masses 
((M 5 oo) = 3, 6, and 9 x 10 14 M Q ) and data types (X- 
ray/simulation, X-ray/SZ, SZ-only) for the A10, P12, 
BoXSZ samples. As another consistency check between 
our data and the results of A10 and P12, we fit each of 
their best-fit gNFW models to our data, allowing only 
the normalization to be a free parameter. Although the 
fit quality is poor in both cases, we find normalizations 
consistent with both results (7.82 compared to 8.40 for 
the A10 fit and 6.54 compared to 6.41 for the P12 fitQ 
This implies that the average total pressure of our sam- 
ple is consistent with the average total pressure found in 
those analyses, further showing the approximate univer- 
sality of cluster pressure profiles and the good agreement 
between SZ and X-ray measurements of those profiles. 
In addition, we note that two other analyses show good 
agreement with the results of A10 (and consequently 
our result s as we ll). These analyses were performed by 
iSun et all (|2011l ) using a set of 43 low redshift groups 

18 Although these single-parameter gNFW fits do not suffer from 
the same degeneracies seen in the multi-parameter fits, we still do 
not include error estimates due to the fact that the fit quality is 
poor. 
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with Chandra data, and by IPlagge et al.l (I2010D using 
SPT SZ data for a set of 15 moderate redshift clusters. 

In addition to fitting gNFW models to our data, we 
also fit a /3-model of the form 



P(X) = 



Pa 



(1 + (X/i? c )2)3/3/2 > 



where Pq is the pressure normalization, R c is 
the core radius, and is the pow er law slope 
(jCavaliere and Fusco-Femianol 119761 119781 ). Fitting the 
full BoXSZ sample, we find best-fit parameters of R c = 
0.11, P = 0.61 and P = 18.9. We find x? cd = L5 for 10 
DOF, which gives a PTE of 0.13 and indicates a some- 
what worse fit compared to the gNFW model. The fit ap- 
pears to be largely driven by the higher S/N data at small 
radii, which explains why the best-fit value of /3 is more 
similar to X-ray derived results f rom fits to the inner re- 
gions of clust ers (/? ~ 2/3, e.g.. iJones fc Formanl [1983 : 
lArnaudl l2009f) than those derived fro m X-ray surface 
brightness profiles at larg e radii (e.g., iVikhlinin et al.l 
ll999l ; lMaugrian et al.ll2008f l or from SZ data extending to 
large radius (J3 ^ 0-85 - 1.05, e.g., lHallman et all [20071 
IPlagge et al1[20To| ). Given the superior fit quality of the 
gNFW model, along with the known shortcomings o f the 
/3-model (e.g.. lMohr et al.lll999t lHallman et aTll2007ft . we 
do not explore the /3-model in any additional detail. 

6. ENSEMBLE PROPERTIES OF THE PRESSURE 
PROFILES 

Jointly fitting distinct subsets of clusters provides some 
information on the differences among cluster pressure 
profiles as a function of radius, but does not directly 
probe the intrinsic scatter among these profiles. Investi- 
gating the scatter instead requires a model for the ensem- 
ble of profiles to be fit to the data from individual clus- 
ters. Here we adopt arguably the simplest such model, 
describing scaled cluster pressure profiles as a Gaussian 
pro cess (for an introduction to G aussian processes, see 
e.g. iRasmussen fc Williamsll2006l ). 

In this approach, the ensemble of profiles is modeled 
by (1) a mean scaled pressure profile as a function of 
scaled radius, P(X)^\ and (2) a covariance function, 
S(Xi, X2), encoding the intrinsic scatter about P(X) 
as a function of radius. Mathematically, the likelihood 
for a single realization of P(X) to have a set of scaled 
pressures {Pi} at scaled radii {Xi} is proportional to 

cxp(-^ T 5" 1 z) 
\S\^ ' 

where Zi = Pi — P{Xi) and Sij = Xj), for any 

choice of radii. Thus, the diagonal covariance terms, 
T,(X,X), dictate the marginal intrinsic scatter among 
profiles at a given radius; while off-diagonal terms, 
E(X 1; X 2 ) with Xi ^ X 2 , determine whether realiza- 
tions of the profile tend to be shifted coherently with re- 
spect to P (positive values) , or tend to cross P (negative 

19 Note that this ensemble mean profile is conceptually different 
from the average parameterized profiles fit in Scction[5] The former 
describes the average profile accounting for the presence of intrinsic 
scatter, i.e. the center of an ensemble of profiles at a given radius, 
while the latter assumes that all cluster profiles are described by a 
single function, with residuals between the model and data entirely 
due to known measurement uncertainties. 



values). In addition to the mean profile and scatter, we 
fit simultaneously for the mass dependence of the pres- 
sure normalization via the parameter ap in Equation [2"P T 1 

The individually deprojected cluster profiles from Sec- 
tion 0] will be used to constrain this model. However, 
those profiles are each constrained at different scaled 
radii, presenting a significant complication to the anal- 
ysis. To simplify the problem, we interpolate the indi- 
vidual profiles to a set of 5 common scaled radii, log- 
arithmically spaced between 0.1i?soo and 2.0i?5oo- The 
radial range probed here is smaller than that covered 
by the sample as a whole, reflecting the fact that data 
from a sufficient number of clusters must exist at each 
radius in order to constrain the scatter. The interpo- 
lation was accomplished as follows. First, we generate 
a multivariate normal draw of the pressures at the de- 
projection radii using the measured mean values and 
measurement error covariance matrices of the individ- 
ual cluster deprojections from Section H) We then in- 
terpolate these pressures to the set of 5 common scaled 
radii using a power-law interpolation (recall that our de- 
projections assumed power-law behavior between the dc- 
projection radii). This process is repeated many times 
in order to constrain the mean pressures and the mea- 
surement error covariance matrices at the common scaled 
radii. Since the individual profiles do not cover the entire 
range 0.1 < X < 2, each cluster provides information at 
only a subset of the final radii. The result of this proce- 
dure is that the mean profile and covariance function can 
be compactly parametrized by 5 pressure values (for P) 
and the independent elements of a 5 x 5 covariance ma- 
trix (for E), corresponding to the common scaled radii, 
where otherwise we would have been forced to assume a 
particular functional form for T,(X\, X 2 ). 

Thus, the complete log-likelihood used to constrain the 
model is 



zJiS+Uj) Zj +]n\S + Uj 



(3) 



where the sum is over clusters, and Uj is the measure- 
ment error covariance matrix for the interpolated, scaled 
pressure profile of the jth cluster. The parameter space 
for this model was explored using MCMC, adopting flat 
priors on all 21 free parameters. Maximum-likelihood 
confidence intervals for each parameter are displayed in 
Table [6l and Figure [9] shows the recovered pressure pro- 
file and fractional scatter as a function of radius. The 
reduced \ 2 °f our data with respect to the best fitting 
model is 1.02, indicating that the Gaussian process de- 
scription of the ensemble of profiles provides a sufficient 
goodness of fit. In particular, while the difference in 
gNFW profile fits of cool-core and disturbed clusters is 
evident at O.l-Rsoo in Figure[6l its modest statistical sig- 
nificance at that radius is reflected in this analysis, where 
a simple Gaussian scatter is seen to be a satisfactory 
model for our dataFI 

20 Including a free power of E(z) has a negligible effect on our 
results in this section. Since the data constrain this evolution term 
very poorly, we did not investigate such an evolution term further. 

21 Wc also attempted to constrain the intrinsic scatter within 
the disturbed and cool-core subsamplcs. Unfortunately, neither 
subsamplc provides enough data at these small radii for us to con- 
strain ensemble models for the individual subsamplcs. 
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Fig. 9. — Left: the dark, inner shaded region shows the 68.3% confidence posterior for the mean pressure profile determined from our 
Gaussian process analysis, while the light, outer region indicates the best-fit marginal intrinsic scatter at each radius (the square root of the 
diagonal elements of the covariance matrix). At all radii, the uncertainty on the mean function is smaller than the corresponding nominal 
intrinsic scatter. The dashed line shows the best gNFW joint fit, which assumes that a single profile describes all clusters, and differs 
slightly from the Gaussian process fit, which includes intrinsic scatter. Right: the fractional intrinsic scatter (diagonal covarianccs scaled 
by the best fitting mean profile) as a function of radius. At all radii, zero intrinsic scatter is excluded at > 95.4% confidence. The intrinsic 
scatter estimated by A10 from X-ray data is shown as the dot-dash line, and is in good agreement with our results at R < i?soo- 

We constrain the mass scaling to be ap = 0.49^; Jg, 
which is shallower than the value of 2/3 predicted by 



self-similar hydrostatic equilibrium scalings Laiserlll986[ 
A10). Under the assumption that the pressure profile 
shape is independent of mass, the integrated SZ signal 
Y, or its X-ray analog Yx, scale with mass according to 
a power law slope ay = ap + 1. Our value of ap there- 
fore implies a Y-M scaling with a power law slope of 
1.49, in good agreement with the Yx~M power law slope 
of 1.48±0.04 found by M10 using an identical method for 
constraining cluster masses from Chandra X-ray dataPl 
These results are inconsistent with those of A10, who 
measured a Yx~M scaling with ay = 1.78 ± 0.06, im- 
plying ap = 0.78. We speculate that this discrepancy 
in values of ap is due to differences in mass estimation 
between A10 and our work (whic h uses mas s es det er- 
mined according to M10). Indeed. iRozo et all (|2012af bT) 
have noted that the difference in Yx~M slopes between 
A10 and M10 is consistent with being due to a system- 
atic disagreement in mass estimates, after accounting for 
instrumental calibration. Comparable disagreements in 
X-ray temperature-mass relation slopes, loosely corre- 
lated with hydrosta tic mass calibrat i on tec hniques, were 
also pointed out bv lMantz fc Allen"! (| 2 1 If ) . In addition, 
as described in Section [4j we emphasize that our results 
for the shape of the average pressure profile arc approxi- 
mately independent of the exact mass scaling, and there- 
fore should not be affected by any of these discrepancies 
in mass determinations. 

Our fit detects the presence of non-zero intrinsic scat- 
ter among profiles (diagonal covariance terms) at a sig- 
nificant level (> 95.4% confidence) at all radii, with the 
fractional intrinsic scatter minimized at radii ~ 0.2i?5nn- 
0.5i?5oo- The off-diagonal covariance terms become con- 
sistent with zero at large radial separations (i.e. large 

22 A separate analysis of the BoXSZ sample, which accounts for 
selection effects and uses directly integrated cylindrical Y values 
wit hin -R2500, also finds a similar slope for the scaling of Y versus 
M JCzakon et al. in p repl ) . 



lA'i — A2I), indicating that our pressure scaling has ef- 
fectively removed the principal mass and rcdshift depen- 
dence of the profiles. The adequacy of our simple scal- 
ing over a wide redshift and mass range further confirms 
that cluster pressure profiles are approximately universal, 
with fractional intrinsic scatter away from the universal 
profile at the ~ 20-40% level. 

In addition, given our measurement uncertainties, the 
intrinsic scatter of our high mass and moderate redshift 
BoXSZ sample is consistent with that found using X-ray 
data within -R500 for the lower mass and lower redshift 
clusters in the REXCE SS sample (A10 ; Figure U) and 
for low redshift groups (jSun et al.ll2011h . supporting the 
proposition that the intrinsic scatter among pressure pro- 
files is relatively independent of mass and rcdshift. Given 
the large redshift range spanned by the BoXSZ sample, 
this comparison to results at lower reds hifts also indic ates 
that the inexact redshift scalings of the iKaiserl (|1986l ) re- 
lations do not introduce a significa nt amount of intrin- 
sic sc atter among pressure profiles (jKravtsov fe Borganil 
l2012h . Furthermore, similar to the findings of P12, our 
results imply that the intrinsic scatter continues to in- 
crease at radii larger than those probed previously by 
X-ray analyses (> -Rsoo)- 

We note that the Gaussian process approach employed 
in this section provides a compact framework to obtain si- 
multaneous constraints on the ensemble mean profile, its 
scaling relations, and its radius-dependent intrinsic scat- 
ter. In the present analysis, we have made some assump- 
tions for computational expediency, namely that the 
measurement errors and intrinsic scatter are described 
well as Gaussian; however, neither of these assumptions 
is necessarily required. Looking forward, such probabilis- 
tic models in general are attractive because they allow 
straightforward and quantitative comparisons of clusters 
from different observed samples or from simulations. For 
example, given such a description of simulated cluster 
profiles, it would be straightforward to test whether ob- 
served profiles are consistent with the simulations using 
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TABLE 6 

Gaussian process fit parameters 



radius (R 50o ) 0.100 



0.212 



0.449 



0.951 



2.013 



mean scaled pressure profile, P(X) 



13.3 



+ 1.8 
1.8 



4.59 



+0.20 



1.47 



+0.08 



intrinsic scatter covariance matrix, S(Xi, X2) 



0.100 
0.212 
0.449 
0.951 
2.013 



27.5 



+30.0 
15.0 



2.75 
1.08 



+ 2.50 
2.00 



0.85 



+0.90 
0.80 



+0.45 
0.45 



0.19 1 



0.13 



+0.08 



0.37 



+0.04 



0.22 



+0.45 
0.35 



0.04 1 



0.037^ 



0.019 1 



0.102] 

0.07; 

0.02! 
0.009;! 
0.0038^ 
0.0047 4 



-0.26 
-0.24 



0.04 
0.018 



Note. — Gaussian process fit parameters describing the mean pressure profile 
of our sample and the intrinsic scatter about this profile. The fit also includes an 
overall mass scaling, constrained to ctp = 0.49^'^; see the text in Section [6] All 
of the diagonal elements of the covariance matrix are different from zero at greater 
than 95.4% significance, indicating that we detect non-zero intrinsic scatter at all 
radii. Most of the off-diagonal elements are consistent with zero, indicating that our 
mass and redshift scaling largely accounts for any evolution in the normalization of 
the pressure profiles. 



simple x 2 statistics. 

7. SUMMARY 

We have examined the pressure profiles determined 
from Bolocam SZ observations of the 45 clusters in the 
BoXSZ sample. This sample spans a large range in 
redshift (0.15 < z < 0.89), with a median redshift of 
z = 0.42. These clusters are also among the most mas- 
sive known, with a median mass of M500 = 9 x 10 14 M Q . 
All of these clusters have Chandra X-ray observations, 
and we have used these X-ray data to determine the mass 
of each cluster. Using these masses, we have scaled each 
SZ pressure profile by the mass-and-rcdshift-dependent 
normalization factor P500 and by the overdensity radius 
-^500- We constrained the average pressure profile of the 
BoXSZ sample using a joint deprojection technique, us- 
ing 13 radial bins approximately logarithmically spaced 
between 0.07i?soo and 3.5i?5oo- The X-ray data were also 
used to classify the disturbed and cool-core subsamples 
of the BoXSZ, and we deprojected the average pressure 
profiles of these two subsamples into 7 radial bins span- 
ning the same radial range (0.07i?5oo to 3.5i?5oo)- 

We fit gNFW models to all three of these average de- 
projected pressure profiles, and we find that this model 
describes our data with a sufficient goodness of fit when 
4/5 of the gNFW parameters are allowed to vary. The 
best-fit average pressure profile of our full sample is de- 
scribed by the parameters (C500, a, P, 7, Po = 1-18, 0.86, 
3.67, 0.67, 4.29). We find a worse, although acceptable, 
fit quality using a /3-model, but we do not explore an 
in-depth analysis of /3- model fits due to the known short- 
commgs of that model (jMohr et al.lll999t lHallman et~aTl 
120071 ). The gNFW fits show consistent pressure pro- 
files regardless of cluster morphological classification out- 
side of ~ 0.15i?5oo given our measurement uncertain- 
ties, but inside that radius the cool-core systems show 
higher pressures than the disturbed systems. Due to 
the large parameter degeneracies in the gNFW model, 



our best-fit parameter values are not in general similar 
to those found in previ ous analyses (iNaeai et al.l 120071 : 
IPlagge et al.ll20T(iL A10,|Suri_eiaIJ|2011] P12). However, 
the actual profile shapes are quite similar, although our 
data provide hints of slightly higher pressures at both 
the smallest and largest radii. This agreement provides 
further evidence that the average cluster pressure profile 
is approximately universal (at least within our measure- 
ment uncertainties, which arc ~ 10-20% within i?5oo), 
especially given the large differences in sample masses 
and redshifts between the BoXSZ and these previous 
analyses. In addition, since many of the previous anal- 
yses relied on X-ray data, rather than SZ data, our re- 
sults show the good ag reement of SZ and X-ray mea- 
surements of the ICM (IPlagge et all 120101: iMelin et all 
20111: \Planck early results xll2011h iKomatsu et alj|2011t 
Bonamente et al.ll2012l P12). 

Finally, we simultaneously fit for the overall mass scal- 
ing, the ensemble mean profile, and the radius-dependent 
intrinsic scatter of the pressure profiles using a Gaussian 
process model. We find that the fractional scatter is min- 
imized at radii between ~ 0.2i?soo and ~ 0.5i?soo at val- 
ues < 20%, with larger scatter at both smaller and larger 
radii. The best-fit mass scaling has a power law slope of 
0.49 (compared to the nominal prediction of 2/3 based 
on self-similar hydrostatic equilibrium models), which is 
nearly identical to the scaling expected from the X-ray 
derived Yx~M scaling determined by M10 using an iden- 
tical X-ray mass determination. Given our measurement 
uncertainties, our intrinsic scatter constraints as a func- 
tion of radius are consistent with previous analyses which 
have largely relied on X-ray meas urements of lowe r mass 
and lower redshift clusters (A10, iSun et aLl[20lTl P12). 
This result provides additional evidence that pressure 
profiles are approximately universal over a wide range 
of masses and redshifts, with intrinsic scatter of ~ 20- 
40% about the universal profile, and that SZ and X-ray 
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measurements of these profiles are consistent with each 
other given current observational precision. 
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